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Abstract. The electronic state around a single vacancy in silicon crystal is investigated by 
^ ' using the Green's function approach. The triply degenerate charge states are found to be widely 

C/3 , extended and account for extremely large elastic softening at low temperature as observed in 

recent ultrasonic experiments. When we include the LS coupling Asi on each Si atom, the 6-fold 
' spin-orbital degeneracy for the state with the valence -1-1 and spin 1/2 splits into Vj doublet 

groundstates and Fs quartet excited states with a reduced excited energy of O(Asi/10). We 
also consider the effect of couplings between electrons and Jahn- Teller phonons in the dangling 
bonds within the second order perturbation and find that the groundstate becomes Vg, quartet 
which is responsible for the magnetic-field suppression of the softening in B-doped silicon. 



1. Introduction 



> 

, Recently, Goto et al. have succeeded in tlie direct observation of the isolated vacancy in high- 

■ purity crystalline silicon (Si), with extremely low vacancy concentration believed less than 10^^ 
cm~^, through the ultrasonic measurements at low temperature [H [2]. They have measured 

^ \ the elastic constants of two type samples, non-doped Si, and boron (B)-doped Si and observed 

■ reciprocal temperature dependence below 20 K down to 20 mK, called 'elastic softening', without 
i the sign of the local distortion [1]. The softening of non-doped Si is independent of the external 

magnetic fields up to 10 T, while that of boron (B)-doped Si is suppressed by the external 
- magnetic fields of the order of 1 T. They have claimed that (1) the softenings are cause by the 

r> \ isolated vacancy groundstates with 3-fold orbital degeneracy coupled to the ultrasound and (2) 

■ non-doped Si has the non-magnetic neutral charge state and B-doped Si has the charge state 
y+ with the valence -|-1 and the spin l/2[lj. Based on the standard theory of vacancyjH SJ [5] 
established in 1980's, the groundstate degeneracies for both and states are resolved with 
symmetry lowering due to the Jahn- Teller distortion in contrast to the newly observed softening. 

To understand such anomalous softenings, new approaches focusing on the quantum many- 
body effect have been performed on the basis of cluster model for the dangling-bond orbitals in 
a Si vacancy in a couple of years [H [TJ |8]. Yamakawa et a/. jH [7] have revealed that the effect of 
couplings between electrons and Jahn- Teller phonons as well as that of Coulomb interaction in 
the dangling bonds play crucial roles for the electronic states of both and states. Matsuura 
and Miyake[8] have studied the effects of spin-orbit interaction and the Coulomb interaction, and 
have found that the Fg quartet groundstate, which is responsible for the magnetic suppression 



of the softening in B-doped Si, is realized in the state by assuming a negative value of the 
spin-orbit coupling on the vacancy state Ave < 0, although the explicit calculation for Ave has 
not been done so far. 

In the cluster model calculations mentioned above, the effect of the spatial extension of the 
vacancy state, which is important to determine the absolute value of the elastic softening, has 
not been considered. In our previous paper [9], we have investigated the electronic state around 
a single vacancy in infinite Si crystal on the basis of the Green's function approach, and found 
that the T2 triplet vacancy states are widely extended up to 20 A and are responsible for the 
extreme enhancement of the Curie constant of the quadrupole susceptibilities resulting in the 
elastic softening at low temperature. The Curie constant for the trigonal mode is considerably 
larger than that for the tetragonal mode as observed in the ultrasonic experimentsfT]. The effect 
of the spin-orbit interaction, however, has not been considered there. 

The purpose of this paper is to elucidate the effect of the spin-orbit interaction on the widely 
extended vacancy states. For this purpose, we extend our previous study[9j for the case with 
the LS coupling Asi on each Si atom, and then we calculate the spin-orbit splitting in the 
vacancy state Ave by using the Green's function approach. We also consider the effect of 
couplings between electrons and Jahn- Teller phonons in the dangling bonds within the second 
order perturbation. 



2. Model and Formulation 

Our model Hamiltonian consists of the tight-binding (tb) Hamiltonian H^^ with LS-coupling on 
each Si site and a vacancy potential H^c, and explicitly given by 



H — Hfh + H^c, 

^tb = E E E tT'4maCjm'. + E E E AsiMr:^cL.c. 
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where c|^^ is a creation operator for an electron at site i and orbital m (= 3s, 3px, 'ipy, Spz, 
spin a, and cj^^ is that for wave vector k and band s (= 1 16). In first term of eq. the tb 
parameters are written by the Slater-Koster parameters and determined so as to fit the tb 
band energies eks to the LDA band energies [12] for the Si crystal, and explicitly given in Ref [9j. 
The second term in eq. ([2]) represents LS-coupling on each Si site with LS-coupling constant Asi 
which is experimentally determined as Asi = 29 meV|13j. Here, the spin-orbit matrix M is 



with 



M = - 
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(4) 



The vacancy potential H^c excludes electrons from the vacancy site by raising the energy levels 
A for the orbitals belong to the vacancy site. For A — > 00, no electron exists at the vacancy 
site and then an effective vacancy state is realized. 

In the absence of the vacancy (A = 0), the Green's function for the perfect crystal is described 

as 



(5) 




positive value of the spin-orbit coupling on 0.001 0.002 

the vacancy state Ave > 0. 



Figure 2. Vacancy levels -Era (^) ^"^^ 
(□) for Asi = 43 meV together with (O ) 
for Asi = as functions of the inverse of total 
number of k-points l/N^- 



where UasO^) is the eigenvector for the energy band e^s given in eq. ^ and orbital m 
with spin a, a = {m,a). In the presence of the vacancy (A ^ 0), the Green's function 
is obtained by solving the Dyson's equations which can be written in the 8x8 matrix 
representation as Gij = G^^ + G^qAGoj, with the vacancy potential matrix (A)^^ = A5af3, 

where (G^j)q/3 = G^^^ is the Green's function for A = given in eq. ([5]) and (Gjj)^^ = G^^ 
is the corresponding Green's function for A 7^ 0. In the limit A 00, Gij — > with i = 
and/or j = 0, and then Gg^- + GgoAGoj ^ 0. By using this relation, Gij is obtained by G as 

3. Results 

We calculate G^"^ in eq. ([5|) by performing the k summation with up to 20 x 20 x 10 = 4000 
mesh points and obtain Gij. When Asi = 0, a remarkable localized level is found in the band 
gap and the total weight of this state is found to be 6 by summing the contribution to this level 
from all sites. Therefore, we find that this level corresponds to the T2 states with the 6-fold 
spin-orbital degeneracy with the energy Et^- This localized level is occupied by two electrons 
in the state and by one electron in the state, respectively. 

By introducing the LS coupling on each Si atom (Asi = 29 meV), the T2 states split into Ty 
doublet groundstates and Fg quartet excited states with a spin-orbit splitting on the vacancy 
states Ave > as shown in Fig(TJ In FiglJl the energy levels for the vacancy states, Et^, -Erg 
and Ey.j measured relative to the valence band top, are plotted as functions of the inverse of 
total number of k-points Then, we obtain a Ni^. 00 extrapolated value of the spin-orbit 

splitting on the vacancy states as Ave = 2.83 meV which is of order of 1/10 of the spin-orbit 
splitting on Si atom Asi = 43 meV|13]. Therefore, the spin-orbit coupling on the vacancy states 
Ave becomes extremely small as compared to Asi but is positive only due to the effect of spatial 
extension of the vacancy state. 

Finally, we consider the effect of couplings between electrons and nonadiabatic Jahn- Teller 
phonons in the dangling bonds by using the Hamiltonian previously studied in Ref [7], which 
is an extention of the early Schliiter's model [5j. Within the second order perturbation theory. 



we find that the groundstate becomes Fg quartet for a reahstic value of the electron-phonon 
coupling constant. The explicit results will be shown in a subsequent paper [H]. 



4. Summary and Discussion 

In summary, we have investigated the electronic state around a single vacancy in infinite Si 
crystal on the basis of the Green's function approach. When the LS coupling Asi on each Si 
atom is taken into account, the T2 states with the 6-fold spin-orbital degeneracy for the 
state split into doublet groundstates and Fg quartet excited states with a reduced splitting 
energy of O(Asi/10). We have also considered the effect of couplings between electrons and 
Jahn- Teller phonons in the dangling bonds within the second order perturbation and find that 
the groundstate becomes Fg quartet which is responsible for the magnetic-field suppression of 
the softening in B-doped Si. 

In the present study, we have assumed that the Tf^-point symmetry is preserved in a single 
vacancy in infinite Si crystal. Even in this case, the lattice relaxation with keeping the T^- 
point symmetry might take place. Therefore, we have also investigated the effect of symmetry 
preserving lattice relaxation on the spin-orbit splitting in the vacancy states Ave by using 
the first-principle calculation with supercell method, and have found that Ave is positive in 
a single vacancy with an unitcell with 63 Si atoms [H]. To be more conclusive, we need further 
investigation to determine the explicit value of Ave with a larger size of supercell. 

The effect of the electron correlation, which is considered to be important to determine the 
many-body groundstate in a Si vacancy as discussed from the cluster model calculations [6l[71|8], 
has not been considered in this paper. The effect can be discussed by including the selfenergy 
corrections on the basis of the present Green's function approach and will be discussed in near 
future. 
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